Observation of a superfluid Hall effect 
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Measurement techniques based upon the Hall effect 1 are invaluable tools in condensed matter 
physics. When an electric current flows perpendicular to a magnetic field, a Hall voltage develops in 
the direction transverse to both the current and the field. In semiconductors, this behaviour is rou- 
tinely used to measure the density and charge of the current carriers (electrons in conduction bands 
or holes in valence bands) - internal properties of the system that are not accessible from measure- 
ments of the conventional resistance. For strongly interacting electron systems, whose behaviour 
can be very different from the free electron gas, the Hall effect's sensitivity to internal properties 
makes it a powerful tool; indeed, the quantum Hall effects [3] are named after the tool by which 
they are most distinctly measured instead of the physics from which the phenomena originate. Here 
we report the first observation of a Hall effect in an ultracold gas of neutral atoms, revealed by 
measuring a Bose-Einstein condensate's transport properties perpendicular to a synthetic magnetic 
field pi]. Our observations in this vortex- free superfluid are in good agreement with hydrodynamic 
predictions [5 , demonstrating that the system's global irrotationality influences this superfluid Hall 
signal. 



Microscopically, the Hall effect results from the Lorentz 
force F = qv x B experienced by particles with charge q 
and velocity v moving in a uniform magnetic field B. In 
the e x -e y plane perpendicular to B = Be z , this force acts 
on a current with density J = J x e x + J y e y to generate 
an electrochemical potential gradient VV = p H J normal 
to J, where the Hall part of the 2D resistivity tensor 

h - ( ° ~ pxy \ 

is antisymmetric. In conventional metals and semicon- 
ductors, the Hall resistivity p xy = B/qn(r) is related to 
the carriers' charge q and density n(r), but not to the 
dissipative resistivity tensor p = p xx t, where 1 is the 
2x2 identity matrix. Typically, experiments measure 
a sample's longitudinal and transverse voltages V xx and 
V xy (Fig. [Ttt) from which the resistivity tensor can be in- 
ferred [6]. Here, we report an analogous transport mea- 
surement of the full resistivity tensor, including the an- 
tisymmetric contributions from the Hall effect, in a flat- 
tened elongated BEC subjected to a synthetic magnetic 
field B*e z (in which only the charge-field product [7] 
g*B* is defined). 

The transport characteristics of systems with many 
particles and sufficiently strong interactions (Coulomb 
repulsion in electron gases and plasmas, or s- wave con- 
tact interactions in BECs) resemble those of classi- 
cal fluids and are described by hydrodynamics. These 
hydrodynamics can describe ultracold Bose [8 and 
Fermi [9] gases, or characterize the collective modes - 
plasmons [10 and magnetoplasmons [11 - of 2D elec- 
tron gases (2DEGs). We show that a BEC in a synthetic 
magnetic field obeys hydrodynamic equations similar to 
those describing a 2DEG in a uniform magnetic field. 

In a simple Drude model [12] description of a 2DEG 



in a magnetic field (Fig. [Tp) , scattering from impurities 
gives rise to resistance while a position-dependent po- 
tential V(r) controls the electron density. The collective 
electron dynamics can be expressed in terms of hydrody- 
namic continuity and Euler equations: 

= d t n(r) + 1 V • [ Pm (r)n(r)] (1) 

= + ( £ T ^ ) Mr) + V \U(r) + 9n(r)] , 

(2) 

where r = (x,y)\ p m = [m/qn(r)]3 is the mechanical mo- 
mentum of particles with mass m, charge g, and density 
ra (r); = qB/m is the cyclotron frequency for magnetic 
field strength B; U(r) = qV(r) is the potential energy; 
g = 27rh 2 /m accounts for the effects of Fermi pressure 
in this noninteracting 2DEG [10 j; r is the momentum 
relaxation time due to scattering from impurities; and 
D/Dt = dt + 7n _1 [p m (r) • V] is the convective deriva- 
tive. The matrix in equation Q is proportional to the 
resistivity tensor [through the factor m/q 2 n(r)] and has 
equal diagonal components from p (conventional resis- 
tance) and antisymmetric off-diagonal components from 
p H (Hall resistance). To compare trapped interacting 
BECs with the 2DEGs described above, the superfluid 
hydrodynamic (SFHD) equations may be extended [5|fT3] 
to include a time-independent artificial magnetic field 
B* = B*e z . The resulting 2D continuity and Euler 
equations are identical to equations and ([2|, com- 
plete with the off-diagonal antisymmetric components of 
the resistivity tensor - the Hall term - proportional to 
the cyclotron frequency Qq = q*B* jm (see Supplemen- 
tary information). The BECs trap provides the external 
potential energy U(r) = m(u 2 x 2 + u 2 y 2 )/2 with trap 
frequencies uj XiV along e XjV . The interaction term gn(r) 
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FIG. 1. Hall effect measurements, a, 2D electron gas 
model, pictured as a Hall bar, depicting the transverse (Hall) 
and longitudinal voltages V xy and V xx . At the microscopic 
level (right), the electrons move in curved paths under the 
influence of the magnetic field, and scatter from random im- 
purities (red crosses) with a rate r _1 (Drude model), b, The 
experiment used a 3D BEC, flattened along e z and elongated 
along e y , trapped at the intersection of two 1064 nm laser 
beams (pink). A nearly uniform synthetic magnetic field 
B* = B*e z (green arrow) was created via the combination 
of two Raman laser beams (red arrows) counterpropagating 
along e x , and a magnetic field gradient along e y , making a 
vector potential A* = (— B*y, 0,0). c, We generated a cur- 
rent in the BEC by modulating the amplitude of the external 
trapping potential along e^. The BEC was alternately re- 
laxed (top) and compressed (bottom). The thick black lines 
schematically indicate the extent of the BEC in the absence 
of any current and the blue ellipses indicate the regions oc- 
cupied by the BEC as the system was modulated. When 
B* — (left), the transport was primarily along e x , while 
when B* ^ (right), the current was deflected along e y : a 
superfluid analogue of the Hall effect. 

accounts for the mean-field energy due to contact inter- 
actions in the BEC. In principle, the resistivity tensor's 
diagonal components are zero due to the BEC's super- 
fluidity; however, atom loss can imitate damping and 
we retain a phenomenological dissipative term propor- 
tional to t -1 . Though the systems with which we work 
are fully three-dimensional, the dynamics of interest are 
in the plane perpendicular to g*B* and are described 
by these 2D equations. For quantitative comparisons to 
data, calculations were performed using the 3D version of 
equations ([!]) and ([2| (see Supplementary information). 

By expressing the SFHD equations in terms of the 
gauge- invariant mechanical momentum p m , which de- 
scribes the actual motion of particles, we depart from 
the conventional expressions in terms of the canonical 
momentum p = p m + <?*A* (where A* is the artificial 



vector potential such that B* = V x A*), which is di- 
rectly related to the gradient of the superfluid order pa- 
rameter's phase. The distinction between p and p m has 
tangible consequences when considering the condition of 
irrotationality V x p = 0. When a synthetic magnetic 
field is present, irrotationality requires equilibrium fluid 
flow, with V x p m = — g*B*. This need not be the case 
at places like vortex cores, where superfluid density is 
zero; however, in our trapped BEC, vortices are ener- 
getically allowed only when f^c is greater than a critical 
value [U [14]. Unless noted, we remained in the irrota- 
tional regime where is below this critical value. 

In this experiment, we studied mass transport in a 
BEC subject to a synthetic magnetic field B* created 
using Raman dressing [4j [7]. While previous experi- 
ments used rotating traps to exploit the equivalence be- 
tween the Coriolis and Lorentz forces [15] , methods that 
use atom-light coupling [16] instead of rotation allow for 
greater geometric freedom and potentially increased field 
strengths [17]. The Raman technique used here enabled 
our experiment's highly anisotropic BECs, whose 4:1 as- 
pect ratio (Fig. [IJd) mimicked a typical bar-like Hall ge- 
ometry and facilitated identification of the Hall response. 

Our experiments began with a nearly-pure 87 Rb BEC 
of about 2 x 10 5 atoms in the \F = = —1) inter- 

nal state [18] confined in a crossed optical dipole trap 
with frequencies (uu x , uo y , uj z ) /2tt « (11.3,42.5,90) Hz 
(Fig. [TJd). Laser beams counterpropagating along e x 
Raman-coupled spin states rriF=i = 0, ±1 of differing 
momentum to create new dressed eigenstates - spin- 
momentum superpositions [4 . The dressed atoms ex- 
perienced a synthetic vector potential [4], which in the 
presence of an appropriate real magnetic field gradient 
gave g*A* = (— q*B*y, 0, 0) and the associated synthetic 
charge-field product q*B* = q*B*e z . 

To measure transport, we generated a mass current 
(a velocity field) in the BEC by perturbatively modu- 
lating the trap along e x (Fig. [1J). We modified the 
potential energy U(x,y,z) = {n x x 2 + n y y 2 + k z z 2 )/2, 
where ^ = m^f is the spring constant along e^, by 
adding a time-dependent drive 5U X = (5k x x 2 /2) sin(cjt) 
with amplitude 5n x and drive frequency uj. To ensure 
linear response, we chose 5n x < 0.13^; to avoid tran- 
sients, we smoothly increased 5k x from zero in 500 ms. 
After driving in steady-state for a time t mo d, we measured 
the atomic density, either in situ or after time-of-flight 
(TOF). We characterized the system's temporal response 
to the drive by repeating this procedure for 64 equally 
spaced values of t mo d, the longest of which encompassed 
eight cycles of modulation. 

When the geometry of the trap deformation is well- 
overlapped with one of the system's low-energy collective 
excitations [19j [20] - eigenmodes of the SFHD equations 
- the response is enhanced as uj approaches that mode's 
resonance frequency. Such collective modes have been 
measured in both Bose and Fermi gases [8, 9 . In this 
work, we predominantly drove the quadrupole-like "J2" 
eigenmode (whose frequency is interaction-dependent) 
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FIG. 2. Time-dependent Hall signal, a-f In situ ab- 
sorption images reveal the atomic density distribution n(x,y) 
at different times during the 83.3 ms modulation cycle: a,d 
t =157.3 ms, b,e t =178.1 ms, c,f t =198.9 ms. For tt c = 
(a-c), the density modulation is primarily along e x , while for 
Qc/27r = 12.4 Hz (d-f) there is a visible skew to the cloud, 
which changes as a function of modulation time. g,h The 
evolution of the second-order moments of n(x,y) [equation 
|3|] during the steady-state modulation of the trapping po- 
tential, without (g) and with (h) synthetic magnetic field. 
Pale curves are guides to the eye. 

corresponding to compression and expansion primarily 
along (an oscillatory (x 2 ) response), with smaller 
out-of-phase contributions, due to interactions, along e y 
and e z ({y 2 ) and (z 2 ) responses). When the synthetic 
magnetic field is non-zero, correlated transport along e x 
and e y appears (an (xy) response) resulting from mixing 
with the scissors mode [21] . This independent corre- 
lated transport is the Hall effect. Had we excited the 
"Xl" dipole (sloshing) mode, we would have observed 
the Hall effect, but would not have probed internal prop- 
erties related to interactions, as described by Kohn's the- 
orem [22]. 

We absorption-imaged the time-dependent column 
density distributions n(x, y) either in situ, to directly de- 
termine the density response, or after TOF, to determine 
the momentum response. For in situ measurements, we 



FIG. 3. Spectra of response to trap modulation. These 
responses, measured after TOF, are dominated by the in situ 
momentum. Second order moment amplitudes A(x 2 ), A{y 2 ), 
and A(xy), as a function of modulation frequency to (a, b, and 
c, respectively). Top row, In- and out-of-phase responses 
(open and closed symbols, respectively), as determined by 
a lock-in analysis (see Supplementary information). Dashed 
(solid) curves show the results of a SFHD fit to the response 
for the in-phase (out-of-phase) components. The fit was per- 
formed simultaneously to all six sets of data. The uncer- 
tainties reflect the standard deviation of the mean. Bottom 
row, Amplitude response (root-sum-square of the above re- 
sults). Solid curves show the calculated response using the 
above fit results. 

imaged the atoms immediately following release from the 
trap (t ff < 1 /is). For TOF measurements, we allowed 
the atoms to expand freely for 36.2 ms before imag- 
ing. During the first 2 ms of TOF, we deloaded [23] 
the Raman-dressed atoms into the \F = 1, rap = 1) bare 
state before removing the Raman dressing fields. 

These density distributions offer snapshots of the dy- 
namically evolving BEC at t mo d- For each t mo d, we 
extract three independent second-order moments of the 
density distribution, (x 2 ), (y 2 ), and (xy): 



{X{X j 



Y.XjXj n(x,y) 
^n(x,y) 



(3) 



where Xi G {x, y} was measured from the centre of the 
density distribution ((x) = (y) = 0). 

Figure [2] shows n(x,y) at three points during the modu- 
lation cycle without (Figs.[2|i-c) and with (Figs.[2]i-f) the 
synthetic magnetic field. While the compression of the 
cloud along the drive direction e x was always the most ev- 
ident signal, a nonzero synthetic magnetic field (Qq 7^ 0) 
tilted the clouds as the velocity field responded to the 
drive current, macroscopically manifesting the BEC's su- 
perfluid Hall resistivity. In Figs.[2]i-f, the Hall signal (xy) 
oscillates about zero and reaches its minimum, zero, and 
maximum, respectively. The temporal responses of the 
second-order moments are plotted in Figs. [2^ (Qc = 0) 
and [2J1 (Qq 7^ 0) f° r oj/2tt = 12 Hz. This modulation 
frequency is above resonance (Fig. [3]): close enough for 
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FIG. 4. Hall effect, a, Measurements of the Hall signal 
(xy) I R x R y (black) for modulation frequency cj/27t = 12 Hz 
at various synthetic magnetic field values, from in situ im- 
ages. Thomas- Fermi radii R x ,R y are from fits to the data. 
Filled circles indicate data in which the system has no vor- 
tices and open circles indicate data for which TOF images 
indicate that vortices were likely present. The top axis indi- 
cates the calculated "magnetic flux" penetrating the BEC [4], 
using &b* /$o = QcmA/h, where A = ixR x R y is the area en- 
compassed by the BEC in the e^-e^ plane, which increases 
with Qc • The vertical uncertainties reflect the standard devi- 
ation of the mean; the horizontal uncertainties reflect only 
the experimental uncertainty in Qc and not the Thomas- 
Fermi radii, b, The synthetic field dependence of the time- 
dependent part of the normalized second order moments 
A(x 2 )/R 2 x (blue circles) and A(y 2 ) / ' R 2 y (red circles). The solid 
curves show the calculated SFHD response, which was simul- 
taneously fit to all three moments. The signals' nonlinearity 
is due to the decreasing collective mode resonance frequency 
with increasing Qc- 



a strong signal, but far enough to spectrally resolve the 
natural mode from the driven mode. While (x 2 ) and 
(y 2 ) were not significantly affected with the introduction 
of the synthetic field, a weak f^c-dependence of the par- 
allel (x 2 ) and perpendicular (y 2 ) moments arises from 
variations in the effective trapping frequencies as field 
strength increases [7J [13] . 

Figure [3] shows the spectrum of the response across 
the low- frequency X2 resonance. We measured the time- 
dependent second-order moments (x 2 ), (y 2 ), and (xy) 
and used a lock-in technique to determine the response 
amplitudes A(x 2 ), A(y 2 ), and A(xy) for each moment 
both in-phase and out-of-phase with the drive. We 
used the TOF method to ensure a strong signal even 



far from resonance. At a moderate synthetic magnetic 
field strength [Qc/27r = 10.5(5) Hz], no vortices entered 
and equations 0, ([2| apply. We simultaneously fit all 
six cj-dependent responses to the solutions of the lin- 
earized SFHD equations (propagated through TOF with- 
out linearization, see Supplementary information) and 
extracted three free parameters: the phenomenological 
damping term r _1 = 4.4(4) s _1 , which is the diagonal 
part of the resistivity matrix in equation ([2]); the effec- 
tive trap frequency, uj x /2tt = 6.07(2) Hz; and the drive 
strength ^k x )k x — 0.074(2). The latter two parameters 
are in reasonable agreement with the directly measured 
quantities u x /2tt = 6.1(2) Hz and 5k x /k x = 0.10(1). 

While the above data are well-characterized by and 
are in good agreement with a SFHD description and its 
TOF propagation, in situ measurements provide a clearer 
demonstration of Hall physics. For different values of 
^c, we measured the BECs' time-dependent density pro- 
files and calculated normalized second-order moments 
(x 2 )/R 2 , (y 2 )/R 2 , and (xy) / R x R y from fits to skewed 
Thomas-Fermi distributions with radii R x and R y (see 
Supplementary information) . As in Fig. [2J the modula- 
tion drive frequency was uj/2tt = 12 Hz. 

Figure [4] shows the amplitude response of the normal- 
ized second order moments along with simultaneous fits 
to the linear-response predictions of SFHD using a sin- 
gle fit parameter, 5n x (with the previously measured uj x 
and t -1 as fixed parameters). For small ^c? we find 
good agreement between our measurements and SFHD. 
The nonlinearity of A(xy) / R x R y at high fields is caused 
by the shift of the X2 resonance away from the fixed 
drive frequency uo as VLq increases [7] [13]. For small 
l^c, where the resonance is nearly constant, the slope 
of A(xy) I R x R y vs is a measure of the Hall resis- 
tivity p xy = mQc/Q 2 n averaged over the cloud's chang- 
ing density. The good agreement between our measure- 
ments and the irrotational SFHD model demonstrates 
the utility of this new Hall technique as a simple, quan- 
titative way to probe microscopic properties, like the su- 
perfluid Hall resistivity, of a quantum degenerate gas. 
Beyond VLq/^ ~ 15 Hz, indications of vortices were evi- 
dent in TOF; therefore, we expect that the predictions of 
vortex-free SFHD no longer apply. In the large-!?* limit 
of "diffused vorticity," [24] in which many vortices have 
entered the BEC and "ordinary" fluid behaviour with 
V x p(r) = g*B* is valid, calculations show that the su- 
perfluid Hall signal approaches zero, but does not become 
negative as in some superconducting materials [25j [26] . 

Inspired by condensed- matter techniques, we have 
demonstrated Hall physics in a BEC and showed that 
it is sensitive to the properties of an interacting, irrota- 
tional BEC. This macroscopic Hall measurement tech- 
nique can be extended to probe the microscopic proper- 
ties of more complicated ultracold configurations, such 
as 2D systems with vortices (arising from thermal fluctu- 
ations [27] or synthetic magnetic fields [4] ) , systems with 
spin-orbit coupling [28], and quantum Hall systems [29 . 



5 



[1] Hall, E. H. On a new action of the magnet on electric 
currents. Amer. J. Math. 2, 287-292 (1879). 

[2] von Klitzing, K., Dorda, G. & Pepper, M. New method 
for high-accuracy determination of the fine-structure con- 
stant based on quantized Hall resistance. Phys. Rev. Lett. 
45, 494-497 (1980). 

[3] Tsui, D. C, Stormer, H. L. & Gossard, A. C. Two- 
dimensional magnetotransport in the extreme quantum 
limit. Phys. Rev. Lett. 48, 1559-1562 (1982). 

[4] Lin, Y. J., Compton, R. L., Jimenez-Garcia, K., Porto, 
J. V. & Spielman, I. B. Synthetic magnetic fields for 
ultracold neutral atoms. Nature 462, 628-632 (2009). 

[5] Recati, A., Zambelli, F. & Stringari, S. Overcritical rota- 
tion of a trapped Bose-Einstein condensate. Phys. Rev. 
Lett. 86, 377-380 (2001). 

[6] Simon, S. H. Comment on "Evidence for an anisotropic 
state of two-dimensional electrons in high Landau levels" . 
Phys. Rev. Lett. 83, 4223-4223 (1999). 

[7] Spielman, I. B. Raman processes and effective gauge 
potentials. Phys. Rev. A 79, 063613 (2009). 

[8] Dalfovo, F., Giorgini, S., Pitaevskii, L. P. & Stringari, S. 
Theory of Bose-Einstein condensation in trapped gases. 
Rev. Mod. Phys. 71, 463-512 (1999). 

[9] Giorgini, S., Pitaevskii, L. P. & Stringari, S. Theory of 
ultracold atomic Fermi gases. Rev. Mod. Phys. 80, 1215- 
1274 (2008). 

[10] Fetter, A. L. Electrodynamics of a layered electron gas. 

I. Single layer. Annals of Phyiscs 81, 367-393 (1973). 
[11] Das Sarma, S. & Quinn, J. J. Collective excitations 

in semiconductor superlattices. Phys. Rev. B 25, 7603 

(1982). 

[12] Ashcroft, N. W. & Mermin, N. D. Solid State Physics 
(Harcourt Brace College Publishers, 1976). 

[13] Cozzini, M., Stringari, S., Berin, V., Rosenbusch, P. & 
Dalibard, J. Scissors mode of a rotating Bose-Einstein 
condensate. Phys. Rev. A 67, 021602(R) (2003). 

[14] Lundh, E., Pethick, C. J. & Smith, H. Zero-temperature 
properties of a trapped Bose-condensed gas: Beyond the 
Thomas-Fermi approximation. Phys. Rev. A 55, 2126- 
2131 (1997). 

[15] Fetter, A. L. Rotating trapped Bose-Einstein conden- 
sates. Rev. Mod. Phys. 81, 647 (2009). 

[16] Dalibard, J., Gerbier, F., Juzeliunas, G. & Ohberg, 
P. Colloquium : Artificial gauge potentials for neutral 
atoms. Rev. Mod. Phys. 83, 1523-1543 (2011). 

[17] Cooper, N. R. Optical flux lattices for ultracold atoms. 
Phys. Rev. Lett. 106, 175301 (2011). 

[18] Lin, Y.-J., Perry, A. R., Compton, R. L., Spielman, I. B. 
& Porto, J. V. Rapid production of 87 Rb Bose-Einstein 
condensates in a combined magnetic and optical poten- 
tial. Phys. Rev. A 79, 063631 (2009). 

[19] Stringari, S. Collective excitations of a trapped Bose- 
condensed gas. Phys. Rev. Lett. 77, 2360 (1996). 

[20] Murray, D. R., Barnett, S. M., Ohberg, P. & Gomila, 
D. Elementary excitations of a Bose-Einstein condensate 
in an effective magnetic field. Phys. Rev. A 76, 053626 
(2007). 

[21] Marago, O. M. et al. Observation of the scissors mode 
and evidence for superfluidity of a trapped Bose-Einstein 
condensed gas. Phys. Rev. Lett. 84, 2056-2059 (2000). 

[22] Kohn, W. Cyclotron resonance and de Haas-van Alphen 



oscillations of an interacting electron gas. Phys. Rev. 
123, 1242 (1961). 

[23] Williams, R. A. et al. Synthetic partial waves in ultracold 
atomic collisions. Science 335, 314-317 (2012). 

[24] Cozzini, M. & Stringari, S. Macroscopic dynamics of 
a Bose-Einstein condensate containing a vortex lattice. 
Phys. Rev. A 67, 041602 (2003). 

[25] Smith, A. W. et al. Pinning and the intrinsic magnetic- 
field dependence of the mixed-state Hall conductivity in 
amorphous MosSi and YBa2Cu3 07_«5. Phys. Rev. B 56, 
R2944-R2947 (1997). 

[26] Huber, S. D. & Lindner, N. H. Topological transitions 
for lattice bosons in a magnetic field. Proc. Natl. Acad. 
Sci. USA (published online, 2011). 

[27] Hadzibabic, Z., Kriiger, P., Cheneau, M., Battelier, B. & 
Dalibard, J. Berezinskii-Kosterlitz-Thouless crossover in 
a trapped atomic gas. Nature 441, 1118-1121 (2006). 

[28] van der Bijl, E. & Duine, R. Anomalous Hall conductiv- 
ity from the dipole mode of spin-orbit-coupled cold-atom 
systems. Phys. Rev. Lett. 107, 195302 (2011). 

[29] S0rensen, A. S., Demler, E. k Lukin, M. D. Fractional 
quantum Hall states of atoms in optical lattices. Phys. 
Rev. Lett. 94, 086803 (2005). 

[30] Pethick, C. J. & Smith, H. Bose-Einstein Condensa- 
tion in Dilute Gases (Cambridge University Press, Cam- 
bridge, 2002). 

[31] Nikuni, T., Zaremba, E. & Griffin, A. Two-fluid dynam- 
ics for a Bose-Einstein condensate out of local equilib- 
rium with the noncondensate. Phys. Rev. Lett. 83, 10-13 
(1999). 

[32] Castin, Y. & Dum, R. Bose-Einstein condensates in time 
dependent traps. Phys. Rev. Lett. 77, 5315-5319 (1996). 

[33] Lin, Y.-J. et al. Bose-Einstein condensate in a uniform 
light-induced vector potential. Phys. Rev. Lett. 102, 
130401 (2009). 

Acknowledgements We appreciate conversations with 
J. V. Porto. This work was partially supported by ONR; 
ARO with funds from both DARPA's OLE program 
and the Atomtronics MURI; and the NSF through the 
Physics Frontier Center at JQI. L. J. L. acknowledges 
support from NSERC; K. J.-G. acknowledges CONA- 
CYT; and M. C. B. acknowledges NIST-ARRA. 

Author Contributions All authors contributed 
to writing of the manuscript. L. J. L. executed the final 
calculations, and led the data taking effort in which the 
remaining authors participated. W. D. P. and I. B. S. 
supervised this work and I. B. S. proposed the initial 
measurements in accordance to preliminary numerical 
and analytic calculations. 

Author Information Correspondence and re- 
quests for materials should be addressed to I. B. S. 
(ian.spielman@nist.gov) . 



6 



SUPPLEMENTARY INFORMATION 

Drude model hydrodynamics 

In a 2DEG, the equation of motion for the electrons is 
given by Euler's equation [TO] [11] 

Vt r m 

+W(r) + - VP(r) = 0, 

n 

where p m (r,t) = m dr/dt is the electrons' mechanical 
momentum, q their charge, r the mean time between 
scattering events, U(r) = qV(r) the external potential 
energy, and P(r) the local pressure. This hydrodynamic 
description involves the momentum distribution p m (r,t) 
of the fluid's velocity field, rather than the momenta of 
individual particles, requiring the use of the convective 
derivative D/Dt = d t +m~ 1 [p m -\7]. For a noninteracting 
single-component 2DEG, the internal energy in an area A 
is /Tint = 7rh 2 N 2 /mA = NE F /2, where N is the number 
of electrons and E F the Fermi energy. From the energy, 
the local pressure P = — (dU\ n t/ 8A)n = 7rh 2 n 2 /m can 
be determined. Because P oc n 2 , the force due to the 
internal Fermi pressure is proportional to the gradient of 
the density, and is (coincidentally) directly analogous to 
the effect of contact interactions in a BEC. 



Superfluid hydrodynamics 

The SFHD equations described in the text were de- 
rived from the time-dependent Gross-Pitaevskii equation 
(GPE) [30] 



+^(r)|^(r,t))|V(r,t) 



with the order parameter (or wavefunction) ip(r,t) and 
inverse mass tensor m _1 (which may in general be 
anisotropic). By writing the order parameter i/j(r,t) = 

^^ e iu(r,t) e -iq* f A*(r,t)-ds/h [n termg q{ the den _ 

sity n(r, t), the phase f(r,£), and the vector potential 
q*A*(r,i) (where the integral is evaluated over the posi- 
tion of the particle s), and neglecting the quantum pres- 
sure term [30] , the GPE can be reformulated as the con- 
tinuity and Euler equations: 

= d t n(r, t) + m- 1 V • [p m (r, t)n(r)] (M.l) 
= <9 tPm (r, t) + V [|p m (r, t) • m" 1 • p m (r, t) + U (r, t) 
+^n(r)]+r- 1 p m (r,t), (M.2) 

where p m (r, t) = p(r, t) — g*A*(r) is the mechanical mo- 
mentum and p(r, t) = ftV^(r, t) is the canonical momen- 
tum. A phenomenological damping constant r is also in- 
cluded, which is attributed to loss of atoms from the con- 
densate due to photon scattering from the Raman beams. 



It can be shown [31] that such loss, through exchange of 
atoms with the thermal cloud, results in effective momen- 
tum damping within the linear approximation (below) 

In this context, the important property of superfluidity 
is its irrotationality, and this condition, Vxp(r,t)=0, 
places a further constraint upon the dynamics. When 
q*B* is sufficiently small [4j [14] and the Landau gauge 
vector potential A*(r) = (—£?*?/, 0,0) is assumed, the 
equilibrium condition (determi ned by se tting the time 
derivatives to zero in equations ( |M.l ) and (M.2)) is found 
when p = \q*B* (l + e )V (#?/). The solutions of equa- 
tions (M.l) and ( M.2| ) in equilibrium yield the time- 
independent in situ density distribution of the BEC, 
which is of the usual Thomas-Fermi form [30 , but with 



effective trapping frequencies uj 2 



|fic(l 



(jj. 



2 _ 



u 2 + |^c(l — £). The anisotropy parameter 
e = (uo 2 —uj 2 )/(uj 2 +c2^) is determined implicitly [13] for 
each Qc The (^-dependence of these effective trapping 
frequencies uo x and uj y shifts the resonance frequencies of 
the effective modes as a function of synthetic magnetic 
field strength q*B*. 



As written, equations (M.l) and (M.2) apply in any 



dimension. Assuming 2D, a diagonal uniform mass ten- 
sor (m = ml) and th e spa tially depen dent 2D interac- 
tion term, equations (M.l) and (M.2) are transformed 



using standard vector identities to equations (1) and 
(2). In order to compare calculations using the SFHD 
equations to our data, we used the three-dimensional 
(3D) form, with the mean-field interaction parameter 
g = 47rh 2 a s /m (where a s is the s-wave scattering length), 
and the anisotropic effective mass tensor 

1/ra* 
1/ra 
1/ra 

where ra* is the effective mass calculated from the mod- 
ified dispersion relation along e^, which arises from Ra- 
man dressing [7]. 

To calculate the dynamics under a time-periodic mod- 
ulation of the trapping potential, we assumed a pertur- 
bative change 5U X of the trapping potential energy U as 
described in the main text, and for convenience used a 
complex time dependence: U (r, t) = U (r) + 6U x (r)e~' lu;t . 
Unde r this appro ximation, the linearization of equations 
flM.ip and flM^D gives 



=d t 5n(r, t) + m" 1 V [n (r)£p(r, t)+ 

(p (r)-q*A*(r))5n(r,t)] (M.3) 
=dtSp(r,t) + V [(p (r) - g*A*(r)) • m" 1 • 5p(r,t) 

+5U(r) + g5n(r, t)] + T^Spfa t), (MA) 

where we modelled the time-dependent density and mo- 
mentum as the equilibrium distributions plus a small 
contribution oscillating at the drive frequency: n(r, t) = 
no (r) + 5n(r, t) and p(r, t) = po (r) + 5p(r, t) . (Note that 
for our experiments 9tA*(r, t) = 0.) In the 3D calcula- 
tions, the density modulation was assumed to be of the 
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form Sn(r,i) = e tUJt (Sn + Sn x x 2 + fe^ 2 



£n 2 



Srixyxy) and the momentum modulation 5p(r,t) = 



2u;t V(^(r) where 5u(r) = (^ + ^^ 2 + ^2/ 2 + < s >^ 2 - 



5v xy xy [which ensures irrotationality: V x 5p(r, t) = 0]. 
When the synthetic field is present, non-zero coefficients 
of the xy terms are required to obtain a solution and 
indicate a Hall response. To extract the temporal depen- 
dence of 5n(r, t) and 5p(r, t), we calculated the responses 
at equally spaced t G [0,2tt/lj], taking the imaginary part 
of [— 5n(r, t)] and [— 5p(r, t)] at each t to obtain the re- 
sponse to the sin(cjt) drive. 

To compare these results with our TOF measurements, 
we propagated the calculated in situ density and momen- 
tum distributions throu gh 36 .2 ms TOF [HI [32] using the 
SFHD equations and M3 with A*, U and r~ x 



set to zero. We projected these numerical solutions for 
the density profile into a two-dimensional plane to mimic 
absorption imaging and analyzed them in the same man- 
ner as the data. 

This theoretical model is fit to the experimentally de- 
termined second-order moments (Fig. 3). As an addi- 
tional check, we compared the experimental results de- 
termined by direct integration with results obtained by 
fitting Thomas- Fermi profiles to the density distributions 
n(x, y) and comparing the fit parameters to those ex- 
pected from the SFHD model. Both methods of analysis 
were in excellent agreement with the model, and each 
other. 



System preparation 

After standard laser cooling and trapping proce- 
dures [18], our experiments began with nearly pure 
87 Rb BECs of approximately 2 x 10 5 atoms in the 
\F = l : rriF = —1) state. The condensate was trapped 
at the intersection of two 1064 nm laser beams, propa- 
gating along e x (with w 500 mW focussed to 1/e 2 radii 
of w y ~ 120 fim and w z ~ 50 /im) and e y 125 mW 
with w x « w z « 160 /im), crossing near their focii. 

We implemented a synthetic vector potential using 
Raman dressing [4j [33]. A (real) uniform magnetic 
field B = B$e y Zeeman split adjacent mp levels by 
E z /h 3.25 MHz. Two Raman beams at A « 790.1 nm 
counterpropagated along (with frequencies 00^1^2) 
and coupled the three mp levels with a strength of 
MIr = 5.6(6)i?L for the TOF measurements and £7r = 
6.3(3)£^l/^ for the in situ measurements, where E\, = 
ft 2 /2mA 2 = h x 3.67 kHz is the recoil energy. The detun- 
ing hS = ^(cjri — 0JR2) — Ez from Raman resonance was 
set to zero at the centre of the trap. To realize a synthetic 
magnetic field, we introduced a spatial dependence to the 
detuning with a (real) magnetic field gradient along e y 
such that B(r) ~ ~Bo + b'ye y ; values of ji^b' /h range from 
to 640 Hz/um in this experiment. The resulting vector 
potential A* = B*y e x gives a synthetic magnetic field 
B* = B*e z , with a cyclotron frequency Qq = g*£?*/ m *- 
The Raman dressing in this geometry altered the disper- 



sion relationship only along e^, giving an effective mass 
m*/m 2.6(1). 



Generating currents 

We modulated the confining potential energy U (r, t) = 
U(r) + 5U x (x) sin(cjt), with amplitude SU x (x) = 5n x x 2 /2 
and frequency uj to drive an ac current along e x . To 
allow the system to reach steady-state, we smoothly 
turned on the modulation amplitude in 500 ms from 
zero to 5k x /k x = 0.074(2) for the TOF images, and to 
5k x /kx — 0.151(3) for the in situ images (a greater mod- 
ulation amplitude was required to obtain an acceptable 
signal in situ). This larger modulation corresponded to 
a change of the in situ size of the condensate of < 20% 
peak-to-peak. We verified that the drive was in the lin- 
ear response regime with measurements at 5n x both less 
and greater than the values used here, and found that 
the response was proportional to 5n x . 



in situ image analysis 

After modulating at full-strength for a time t mo d, we 
prepared the system for imaging: in the final 0.5 ms dur- 
ing which the BEC remained trapped, the magnetic field 
gradient was removed and the magnetic bias field was 
aligned along e z for imaging. This process, rapid com- 
pared to the periods associated with trap frequencies, did 
not noticeably affect n(x,y). Upon removing the optical 
trapping potential, a weak 100 us optical pumping pulse, 
detuned « 140 MHz from the F = 1 to F' = 2 tran- 
sition, transferred about 2 % of the atoms (randomly 
sampled from the distribution) to the F = 2 hyperfine 
manifold. These atoms were resonantly imaged in 40 
us on the \F = 2,m F = 2) \F' = 3, m' F 3) cycling 
transition. The peak optical density was generally about 
1.5. 

We obtained the second-order moments described in 
the text (equation (3)) from the 2D density profiles 
n(x,y), shown in Fig. 2. We checked that the ther- 
mal (non-BEC) component of the cloud did not signif- 
icantly contribute to the optical density signal in these 
images by using fits to the appropriate 2D atomic den- 
sity distribution of a mixed (mean-field BEC + ther- 
mal) trapped cloud. Since our results were not affected 
(within experimental uncertainty) by including the ther- 
mal component, we opted to ignore it in favour of a sim- 
pler analysis. We determined the region over which we 
calculated the second-order moments as follows: we fit 
n(x,y) ~ [l — (x/R x ) 2 — (y/R y ) 2 } 3 / 2 , with a 2D gaussian, 
and analyzed those points within the area for which the 
gaussian's value was > 1 % of its peak. Repeating this 
analysis with a 5 % limit gave similar results, confirming 
that the 1 % limit effectively captures the entire BEC. 

For the in situ data presented in Fig. 4, the absorption 
of the optical pumping pulse as it propagated through the 
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cloud caused a small asymmetry to the observed F = 2 
density profiles. To "smooth over" this asymmetry, we 
assume a 3D Thomas-Fermi density profile of the form: 



n(x,y,z,t) = 



1 2 



R x (t) 



y 



Rzit) 



Ry(t)_ 

^xy 

(t)xy^, (M.5) 



where /i is the chemical potential; R x , y , z (t) = Rx,y,z(fy + 
5Rx,y,z(t) is the time-dependent Thomas-Fermi radius of 
the cloud along e XiViZ with equilibrium radius R x , y , z (0); 
and s xy (t) is a skew parameter characterizing the Hall 
response, whose equilibrium value is zero. Fits of equa- 
tion (M.5), integrated over z, to the observed density 



give radii R x , R y and skew parameter s xy . E xact inte- 
gration over the density profile (equation ( M.5| )) gives the 
second-order moments in terms of the Thomas- Fermi pa- 
rameters: (x 2 ) = 4:R%/ 28 — 7 (s xy R x R y ) 2 and (xy) = 

^^xyRxRy I 28 7 (s xy R x R y ) 

Lock-in analysis is performed on the normalized 
second-order moments to find the amplitude response, as 
plotted. We calculated the "in-phase" response by multi- 
plying the time-dependent signals (x 2 ), (y 2 ) and (xy) by 
sin(cjt) and extracting the average of the product to give 
A(x 2 ), A(y 2 ) and A(xy). Similarly, we calculated the 



"out-of-phase" response by multiplying by cos(cjt). The 
displayed uncertainties in these values represent the stan- 
dard deviation of these means. We calculated the ampli- 
tude response of each moment as the root-sum-square of 
the in- and out-of-phase responses. 



A. TOF image analysis 

After modulating at full-strength for a time £ mo d> we 
suddenly removed the optical trap (in < 1 /xs). During 
the following 2 ms of TOF, we adiabatically swept the 
Raman detuning S from to —77 adiabatically trans- 
ferring (deloading) the Raman-dressed spin-momentum 
superposition state to \F = l^mp = +1). During the 
sweep, the subsequent turn-off of the Raman beams, and 
the remainder of TOF, the condensate underwent free 
expansion, initially driven by the mean field. 

We absorption-imaged the expanded atomic distribu- 
tion along e z after 36.2 ms TOF. From the resulting two- 
dimensional density profile n(x,y), we determined the 
second-order moments described by equation (3). The 
discrete sums were analyzed over a region defined as de- 
scribed above. To account only for the dynamics in the 
condensate, we removed the (30 % to 50 %) thermal back- 
ground from the images by means of a gaussian fit to the 
region outside the condensate. We extracted the lock-in 
response as above. 



